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^ . Abstract 

0^ 

The concept of determinant for a linear operator in an infinite-dimensional space is 
addressed, by using the derivative of the operator's zeta-function (following Ray and Singer) 
^ • and, eventually, through its zeta-function trace. A little play with operators as simple as 



±/ (J being the identity operator) and variations thereof, shows that the presence of a non- 
^ ■ commutative anomaly {i.e., the fact that det (AB) ^ det A det -B), is unavoidable, even for 

H ; 

5r ! commuting and, remarkably, also for almost constant operators. In the case of Dirac-type 

operators, similarly basic arguments lead to the conclusion — contradicting common lore — 
that in spite of being det(^ + ivti) = detiff) — im) (as follows from the symmetry condition of 



the ^-spectrum), it turns out that these determinants may not be equal to J detiff)'^ + m?), 
simply because dei[{Ip + im)(fp — im)\ ^ dei(Jp + im) det(ff> — im). A proof of this fact is given, 
by way of a very simple example, using operators with an harmonic-oscillator spectrum and 
fulfilling the symmetry condition. This anomaly can be physically relevant if, in addition 
to a mass term (or instead of it), a chemical potential contribution is added to the Dirac 
operator. 
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1. Introduction 

Many fundamental calculations of Quantum Field Theory reduce, in essence, to the 
computation of the determinant of some operator. One could even venture to say that, at 
one-loop order, any such theory reduces to a theory of determinants. The operators involved 
are 'differential' ones, as the normal physicist would say. In fact, properly speaking, they 
are pseudodifferential operators (\I'DO), that is, in loose terms 'some analytic functions of 



differential operators' (such as ^/l + D or log(l + D), but not! \ogD). This is explained in 
detail in Refs. [|I|]-|[. 

Important as the concept of determinant of a differential or \E'DO may be for theoretical 
physicists (in view of what has just been said), it is surprising that this seems not to be 
a subject of study among function analysts or mathematicians in general. This statement 
must be qualified: I am specifically refering to determinants that involve in its definition 
some kind of regularization, very much related to operators that are not traceclass. This 
piece of calculus — always involving regularization — falls outside the scope of the standard 
disciplines and even many physically oriented mathematicians know little or nothing about it. 
In a sense, the subject has many things in common with that of divergent series but has not 
been so groundly investigated and lacks any reference comparable to the very beautiful book 
of Hardy |Q. Actually, from this general viewpoint, the question of regularizing infinite 
determinants was already addressed by Weierstrass in a way that, although it has been 
pursued by some theoretical physicists with success, is not without problems — as a general 
method — since it ordinarily leads to non-local contributions that cannot be given a physical 
meaning in QFT. We should mention, for completion, that there are, since long ago, well 
stablished theories of determinants for degenerate operators, for traceclass operators in the 
Hilbert space, Fredholm operators, etc. [^ but, again, these definitions of determinant do 
not fulfill all the needs mentioned above which arise in QFT. 

Any high school student knows what a determinant is, in simple words, or at least how 
to calculate the determinant of a 3 x 3 matrix (and some of them, even that of a 4 x 4 one). 
But many one prominent mathematician will answer the question: What is your favourite 
definition of determinant of a differential operator? with: / don't have any, or: These 
operators don't have determinants! An even more 'simple' question I dare to ask the reader 
(which she/he may choose to ask to some other colleague on its turn) is the following: What 
is the value of the determinant of minus the identity operator in an infinite dimensional 
space? Followed by: And that of the determinant YlneN{~^)"' "^ ^^ ^^ actually equal to the 



product of the separate determinants of the plus Is and of the minus Is? 

In this contribution I will point out to specific situations, some of them having become 
common lore already and other that have appeared recently in the literature, concerning the 
concept of determinant in QFT, and I will try to give 'reasonable' answers to questions such 
as the last ones. 

2. Infinite series and (almost) trivial determinants 

The mathematical theory of divergent series has been very fruitful in taming the infinites 
that have appeared in QFT, from the very begining of its conception. Its role is very 
essential, at least in the first stage of the regularization/renormalization procedure. Euler 
and Borel summation methods, and analytic continuation techniques are there commonly 
used. But some difficulties exist that are inherent to the theory of divergent series (see, 
for instance, @]). One of them is the well known fact that, sometimes, by using different 
schemes, different results are obtained. In a well posed physical situation, the 'right' one 
can then only be choosen after experimental validation. Another problem is to understand, 
in physical terms, what you are doing, while performing say an analytic continuation from 
one region of the complex plane to another [^]. This has prevented e.g. the zeta function 
regularization procedure from getting general aceptance among common physicists. 

The situation concerning infinite determinants is even worse, in a sense. There is no book 
on the subject to be compared, for instance, with the above mentioned one by Hardy and we 
see every day that dubious manipulations are being performed at the level of the eigenvalues, 
that are then translated to the determinant itself and elevated sometimes to the cathegory 
of standard results — when not of lore theorems. The first problem is the definition of the 
determinant itself. Let me quote in this respect from a recent paper by E. Witten [0: The 
determinant of the Dirac operator is defined roughly as 

detV = l[Xi, (1) 

i 

where the infinite product is regularized with (for example) zeta function or Pauli-Villars 
regularization. The zeta function definition of the determinant 

det^P = exp [-Cv'm , (2) 

is maybe the one that has more firm mathematical grounds ||^. In spite of starting from 
the identity: log det = tr log, it is known to develop the so called multiplicative anomaly: 



the determinant of the product of two operators is not equal, in general, to the product of 
the determinants (even if the operators commute!). This happens already with very simple 
operators (as two one-dimensional harmonic oscillators only differing in a constant term, 
Laplacians plus different mass terms, etc.). It may look incredible, at first sight, from the 
tr log property and the additivity of the trace, but we must just take into account that the 
zeta trace is no ordinary trace (for it involves regularization), namely: 

tT^V = U-l), (3) 

so that tr^(A + B) ^ tr^A + tr^i?, in general. Not to understand this has originated a con- 
siderable amount of errors in the specialized literature — falsely attributed to missfunctions 
of the rigorous and elegant zeta function method! 

As an example, consider the following commuting linear operators in an infinite-dimensional 
space, given in diagonal form by: 

Oi = diag (1, 2, 3, 4, . . .), 02 = diag (1, 1, 1,1,...)= /, (4) 

and their sum 

Oi + 02 = diag (2,3,4,5,...). (5) 

The corresponding (^-traces are easily obtained: 

trcOi = Ci?,(-1) = -^, tr^02 = Cr{0) = -^, 

13 

trc(0i + 02) = U-l)-l = --, (6) 

the last trace having been calculated according to the rules of infinite series summation (see 
e.g., Hardy [Q). We observe that 

trc(Oi + O2) - tr^Oi - tr^Os = "^ 7^ 0. (7) 

If this happens in such simple situation, involving the identity operator, one can easily 
imagine that any precaution one can take in manipulating infinite sums might turn out 
to be insufficient. Moreover, since the multiplicative anomaly — as has been pointed out 
before — originates precisely in the failure of this addition property for the regularized trace, 
we can already guess that it also can show up in very simple situations, as will now be 
proven, in fact. The appearance of the multiplicative anomaly prevents, in particular, naive 



manipulations with the eigenvalues in the determinant, as reorderings and splittings, what 
a number of physicists seem not to be aware of. 

For warming up, let us calculate some simple determinants with the zeta function method. 
To start with, take 

oo 

Ai = n n. (8) 



n=l 



We have 



oo -1 

Ci(s) = E ^"' = Ci?(s), Ci?'(0) = -- log(27r), (9) 



SO that 



Ai = exp 
a nice result. In the same way, we obtain 



^log(27r) 



27r, (10) 



A,= l[n-' = ^= (11) 

„=i V27r 

— from C2(s) = Cr{—s) — as should be expected. 

Let us now consider the apparently more simple case: 

oo 

A3 = n A. (12) 

n=l 

This poses a problem to the zeta function method, which must be modified somehow to cope 
with such situation. In fact, the corresponding zeta function, 

C3(=^) = E A-^ (13) 

n=l 

has no abscissa of convergence in the complex plane (since the sequence of eigenvalues is 
neither increasing nor decreasing. This can be, however, naturally solved as follows: by 
taking logarithms (what is inherent with the definition of the zeta function method) and 
using again the rules for infinite series, it is plain that the result is 

oo -1 

logA3 = logA El = -:^ log A, (14) 

where the factor in front of log A may be interpreted as the 'zeta measure' of the set of 
positive natural numbers, thus 

oo 

A3 = nA=A"^- (15) 

Jl=l 



This leads, in particular, to the following results: 

oo 

A4 = n 1 = 1"^ = ±1' (16) 

71=1 

for the determinant of the identity operator, J, and 

00 
A,= lli-l) = {-l)-'^=T^, (17) 

n=l 

for the determinant of the operator — /. As it seems clear that the determinant of the 
identity operator should be 1, this tells us (by choice everywhere of the same determination 
of the logarithm in the complex plane) that the determinant of — / is —i and, that of XI, 



the inverse of the corresponding square root of A in (|l^). Notice that, in this way, we 

are starting to build up a set of consistency rules that are reminiscent, in some manner, 

of the corresponding rules for infinite series [^]. More than this, by use of the logarithm, 

all the ordinary rules for infinite series are appliable to the series of logs of eigenvalues, in 

particular, the ones concerning multiplication by a common factor (used before already), or 

of splitting out a finite number of first terms from the series (that is to say, a finite number 

of first factors from the determinant). However, the splitting of an infinite number of terms 

— or of the whole series into two — is not allowed in general. That is, again, the lesson 

we have learned from the existence of the multiplicative anomaly of the determinant when 

evaluated by the zeta function procedure. An additional comment is in order: in dealing 

with infinite series we always take logarithms, and this introduces an ambiguity in the zeta 

function definition of the determinant. This fact is well known [^ |rU| and is common to 

other regularization methods (as Pauli-Villars'), under different disguises. It can duely be 

taken care of by sticking to one and the same determination of the logarithm during the 

whole calculation. 

An apparent problem — or virtue perhaps? — of the zeta function definition of infinite 

determinants is posed by the following example. It turns out that the determinants 

00 00 

Ae = n (-1)'"-"', A, = n {-ir^\ ■ ■ . (18) 

n=l n=l 

and A5 are all different. This originates in 

C6(-l) = 2CM-l,l/2)-l = -i?2(l/2)-l = -^, 

23 

C7(-l) = 4Ch(-1,1/4)-1 = -252(1/4)-! = --, ... (19) 
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and can be interpreted as due to the change of the 'zeta measure' of the number of factors in 
the product leading to the determinant. We could try to avoid this problem by sticking always 
to the most simple characterization of the eigenvalues series (in this case — 1,— 1,— 1,..., 
any reference to n being superfluous). Things are, regretfully, not that simple. Consider the 

determinants A4, A5, and 

00 

A8 = n(-ir' (20) 

n=l 

and try to make compatible the apparently obvious fact that: 



As = VA4A5. (21) 



The determinant Ag can be obtained in three different ways, that yield the same result. 

(i) We have, through the corresponding zeta function, 

00 -1 

a{s) = log(-l 



^ +i- + 0(s) 



7r2s2 ■ 12 
A8 = exp[-C8'(0)|,J = (-l)-^/^^. (22) 

(ii) On the other hand, taking logs as before, from the zeta function measure of the set of 
exponents, we get 



logA8 = log(-l) ^r2 = log(-l)CK(-l), 

n=l 

As = (-1)^«(-^) = (-1)-^/^^ (23) 



(iii) Also, we may instead choose to take derivatives, term by term, in the first expression 
for the zeta function 



C8(^) = E(-i) 



ns 

n=l 

00 00 



C8'(s) = -E^(-l)""^log(-l), C8'(0) = -log(-l)5:n= ^°^^ ^^ 
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n=l n=l -^^ 

As = exp[-Cs'm = {-1)-'^''. (24) 

Remarkably enough, in all three cases we obtain the same result for this determinant. Let 
us now try to fulfill the factorization condition (pT]). We have: 



00 00 



A4 = n 1 = 1"'^' = ±i> A5 = n (-1) = (-1)"'/' = Ti, 



n=l n=l 

00 



As = n (-1)" = (-1)"'/'' = V^4A5 = 1-'/' i-iy'/'. (25) 



n=l 



The only way to fulfill this property ( pi]) is to accept that: 

oo 

ni = -l !! (26) 

n=l 

On the contrary, if we insists (as almost everybody would agree on) that the determinant of 
the identity is 1, then we must give up the property that the determinant of the alternating 
series of eigenvalues 1, —1, 1, —1, ... is equal to the subdeterminant product of the Is, times 
the subdeterminant product of the —Is. This is the most simple reflection one could ever 
have expected to obtain of the multiplicative anomaly of the determinant! 

If we choose to preserve, at any price, the multiplication property of the determinant and 
give sense to the strange fact that det / = — 1, then we do attain compatibility in Eqs. (^Sf ) 
by setting: 



ni = e''^, [](-!) =e'"/^ ^(-l)" = e3'"/^ n(-l)^" = e='*", ... (27) 



Ti=l n=l n=l n=l 



All these are compatible zeta function definitions of the determinant (they can be fixed 
from acceptable roots of 1 or —1 as given by the zeta function exponents) satisfying the 
multiplication rule. However, it is easy to see that this process cannot go for ever (and thus 
eliminate the anomaly): the following dets cannot possibly fulfill the multiplicative property: 

oo oo oo 

n (-1)'"^' = (-1)-^^/^' ^ n (-1)'" n (-i)- (28) 

n=l n=l n=l 

But maybe this is asking too much, as has been observed before, in which case we are still 
left with the compatible (albeit really weird) choice (P?]). 

Anyhow, it is easy to check that we do not find problems in factorizations like 

oo oo oo 

n(-^) = n(-i) n^' (29) 

n=l n=l n=l 



smce 



C(^) = i-irXnis), C'(0) = -I log(-l) - I log(2vr) 

00 00 00 



n i-n) = e-^'(°) = (-1)-^/^ v^ = n (-1) n n, m 



n=l n=l n=l 

or in the more involved one 



00 



n(-ir^=n(-ir n^^ (31) 



71=1 n=l n=l 



because 



C(^) = E(-l)-"^n-^ = $((-l)-^ s), C'(0) = - + $'(0, 1) = Y^ - 2 iog(2^)' 



oo 



n (-i)"n = e-^'(°) = (-1)-^^^ v^ = n (-1)" n ^^ (32) 

n=l n=l n=l 

being $ the polylogarithin function. The factorization of the determinants holds here again, 
in the zeta function prescription, and this fact does not seem to be that immediate, in view 
of the last calculation as compared with what we had before. 

All the above considerations may sound rather trivial, but actually they are not, and 
should be carefully taken into account before proceding with the sort of manipulations of 
the eigenvalues and splittings of determinants that pervade the specialized literature. 

3. The multiplicative anomaly for Dirac type operators 

Consider the ordinary Dirac equation for a massive spinor 

{Ip + im)^ = 0. (33) 

Usually, the determinant of the Dirac operator acting in this equation is obtained by using 
the following argument (see, e.g., 0]): 

r 1 1/2 

det(^ + im) = det(^ - im) = [det(ff)^ + m^)] . (34) 

This comes about from the fact that the spectrum of the Dirac massless operator Ip has 
the following property: if A belongs to the spectrum, then so does —A, that is immediately 
obtained by use of the 75 operator. Then, it turns out that the first det in Eq. (|3^ is a 
product of pairs of the form: 

(A + im){—X + im) = (—A — im){\ — im) = (A — im){—X — im,), (35) 



the last being the pairs appearing in the second det of Eq. ([3^). This is an algebraic 
argument, but there is also the corresponding geometric one, trivial after representing the 
spectral points in the complex plane. Thus, the first equality in (|3^ is proven, and the 
second seems obvious. 

However, due to the existence of the multiplicative anomaly for infinite determinants 
[namely, the fact that, in general, det{AB) 7^ (det A) (detB)], all these formulas, obtained 
by 'simple' manipulation of the eigenvalues, must be set under suspicion and are in need of a 
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rigorous check. Concerning Eq. (|3^, the second equahty cannot be taken for granted, since 
it may turn out that 



det[(^ + im){}P - im)] ^ detiff) + im) det(^ - im). 



(36) 



We will show below that this is indeed the case, in a very simple, parallel example — 
completely under control — that uses as operator some square root of the harmonic oscillator. 
Indeed, consider the square root of the harmonic oscillator obtained by Delbourgo in 
Ref. [jl2|. This example has potentially some interesting physical applications, for it is well 
known that a fermion in an external constant electromagnetic field has a similar spectrum 
(Landau spectrum). Exactly in the same way as when going from the Klein-Gordon to 
the Dirac equation and paying the same price of doubling the number of components {e.g., 
introducing spin), Delbourgo has constructed a model for which there exists a square root 
of its Hamiltonian, which is very close to the one for the harmonic oscillator. It is in fact 
different from the Dirac oscillator introduced by several other authors, corresponding to the 
minimal substitution p ^ p — iaf. The main difference lies in the introduction now of the 
parity operator, Q. Whereas creation and destruction operators for the harmonic oscillator, 
a^ = P ± iX, are non-hermitian, the combinations V^ = P ± iQX are hermitian and 



H' 



{V 



±\2 



P^ + X^TQ = 2HoscTQ. 



(37) 



Notice that the parity term commutes with -ffosc- Doubling the components {(Ti are the 
Pauli matrices) 



d 

P -^ -tai--, X -^ aix, Q -^ (72, 
ox 



the operators V"^ are represented by 



d 



V^ -^ -iai-- ± asx. 
ox 



(38) 



(39) 



In the sequel, we will only consider the operator V = P+. It has for eigenfunctions and 
eigenvalues, respectively. 



i^ni^) 



-le 



-^x'2/2 



/7„_i(x)±/7„(a;)/V2n 



2"+i [n - l)\^ \ [Hn^i{x) T lln{x)H2n 



Xn = ±v 2ra, n > 1, 



ipoix) 



-2;2/2 / I 



An = 0, 



2v^ 



(40) 
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where the Hn{x) are Hermite pohnomials. 

The two operators we shall consider for the calculation of the anomaly are A = V + V 
and B = V — V, V being a real, constant potential with \V\ < y^, that goes multiplied with 
the identity matrix in the two (spinorial) dimensions (omitted here) . Obviously, this V is to 
be identified with the mass m, to make contact with the case of the Dirac equation that we 
had at the begining. 

Notice that V + V and V — V are hermitian, commuting operators. The multiplicative 
anomaly is defined as 

a{A,B) = \ogdet{AB) - logdet A - logdetE = Ca'(O) + Cb'(O) - Cab'(O). (41) 

The zeta function for the operator T> reads 

oo _ 

C^(^) = E V = Eil + (-1)"^] (v^)"' = [1 + (-l)-l2-/\«(s/2), (42) 



n=l 



Cr{s) being the usual Riemann zeta function, which has a simple pole at s = 1. The zeta 
function for the operators P ± \^ is (see, for instance, |[l3|l ) 



CvMs) = E(A^ ± y)" = Us) + E (^^)"^(^ + ^) ci.(n + .). (43) 



Finally, the zeta function for the operator {T> + V){V — V) =V^ — V"^ is given by: 

^"r(n- 

n\ r(s 



Cv^-v<s) = U2s) + Y: ^,IJ, ^ Cv{2n + 2s). (44) 

n=l 



We see that the anomaly in this simple case can be obtained in terms of the derivatives of 
the Riemann zeta function. The important fact is that it turns out to be non-zero: 

a{V + V,V-V) = 2V'^. (45) 

We should point out that this result can be obtained from the Wodzicki formula for the 



anomaly, even if we are working in a non-compact manifold |T3[ . We thus have a very simple 
example of the presence of a non-trivial anomaly for operators of degree one in a space of 
dimension one (spinorial, however). 

One can argue that a mass term will be absorbed by renormalization and will finally yield 
no physical contribution. Notice, however, that the situation is much more general than the 
specific case considered here, which, however, even in its simplicity already accounts for any 
kind of terms not depending on the space-time coordinates. Thus V can represent — aside 
from a mass term im — a constant magnetic field, a finite temperature term, or a chemical 
potential. Anomalous contributions of this kind cannot be absorbed by renormalization (see. 



e.g., [jl4|) and can acquire a direct physical meaning |jT3|, [15 
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